Barkhausen Noise and Critical Scaling in the Demagnetization Curve 
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The demagnetization curve, or initial magnetization curve, is studied by examining the embed- 
ded Barkhausen noise using the non-equilibrium, zero temperature random-field Ising model. The 
demagnetization curve is found to reflect the critical point seen as the system's disorder is changed. 
Critical scaling is found for avalanche sizes and the size and number of spanning avalanches. The 
critical exponents are derived from those related to the saturation loop and subloops. Finally, the 
behavior in the presence of long range demagnetizing fields is discussed. Results are presented for 
simulations of up to one million spins. 

PACS numbers: 75.60.Ej, 64.60.Ht, 75.60.Ch 
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During the last decade, much progress has been made 
in the study of universality in hysteresis and avalanche 
type noise in both experiments and theory. Hysteresis in 
ferromagnets and associated Barkhausen noise have been 
studied in great detail as a conveniently accessible exam- 
ple system. The zero temperature non-equilibrium ran- 
dom field Ising model (RFIM) and recent variants have 
proven very successful in predicting universal scaling be- 
havior in such systems jl], ^, [|. In particular, an under- 
lying disorder induced non-equilibrium phase transition 
has been found, which may be the cause of the broad 
range of power law scaling observed in many systems 
with avalanche-type noise. It is expected to be relevant to 
Barkhausen noise in disordered (hard) magnets [Q. For 
soft magnets, where the long range demagnetizing holds 
are important, the system appears to naturally operate 
at the single domain wall depinning point ||, || , thus 
being self organized critical |jj . 

So far, most of the theoretical studies of this scal- 
ing behavior have focused on the properties of the sat- 
uration hysteresis loop (except for some recent studies 
§, |, H P). In this letter we report history in- 
duced scaling behavior for subloops and in particular the 
demagnetization curve. Ferromagnetic materials with a 
remnant magnetization can be demagnetized by apply- 
ing an oscillating magnetic held with amplitude slowly 
decreasing from a large initial value to zero. Sometimes 
the final state is also termed AC-demagnetized state. The 
oscillating external held with decreasing amplitude takes 
the system through concentric subloops (Fig. |l|). The 
line connecting the tips of the subloops is known as the 
normal or initial magnetization curve, or demagnetiza- 
tion curve. It is to be distinguished from the virgin curve 
which is obtained by thermal demagnetization We 
used the RFIM to simulate the demagnetization curve for 
different amounts of disorder, especially near the critical 
disorder mentioned above. It turns out that this curve 
reflects much of the scaling behavior found for the satu- 
ration hysteresis loop. This may be surprising, since the 
meta-stable states encountered in the demagnetization 
curve are completely different from those of the satura- 
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FIG. 1: The saturation hysteresis loop (dashed line) with its 
corresponding demagnetization curve (solid line) at a disorder 
of R — 2.9J in a 100 3 system. A few of the many subloops 
performed are shown as dotted lines. 



tion loop. In the following we briefly describe the model, 
then derive the necessary exponent relations, discuss the 
numerical results, and their relevance to experiments. 

In the RFIM the magnetic domains of a real magnet 
are modeled by spins = ±1 arranged on a hyper- 
cubic lattice. Nearest neighbors interact ferromagneti- 
cally through a positive exchange interaction J. Variants 
of the model for soft magnets include additional long- 
range interactions due to the demagnetizing field and the 
dipole-dipole interactions. All spins are coupled to a ho- 
mogeneous external magnetic field H(t) that is varied 
adiabatically slowly. To model disorder in the system, a 
quenched random magnetic field hi is added at each site, 
that is taken from a Gaussian probability distribution 
P(hi) = exp(—h^/2R 2 )/^/2TrR, whose width R is a mea- 
sure of the amount of disorder in the system. We write 
the effective Hamiltonian for a system with TV spins 0] : 

H = -JJ2 S ^ - + h i~ JinfM) Si , (1) 

(*i> « 

where Ji n f is an infinite range demagnetizing field , M = 
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J2iSi/N is the magnetization of the system, and (ij) 
stands for nearest neighbor pairs of spins. Initially (at 
time t = — oo), H(—oo) = — oo and all the spins are 
pointing down. Each spin is always aligned with its local 
effective field hf f = JY,(ij) Sj + H(t) + h t -J mf M. The 
external field H (t) is slowly increased from -co until the 
local field, h e ^J of any spin Si changes sign, causing the 
spin to flip Q. Each spin flip can trigger coupled spins 
to flip as well, thus leading to an avalanche of spin flips. 
During such an avalanche the external magnetic field is 
kept fixed (adiabatic case). It is only changed in between 
successive avalanches, to trigger the next spin flip that 
seeds a new avalanche. The simulations are based on 
code both available on the web |l4| and that obtained via 
correspondence ]l5[ which has been modified to allow for 
subloops in the history. In the following we set Jj n / = 0. 

A similar process to the AC demagnetization of the 
laboratory is used in simulating the demagnetization 
curve. For a particular realization of random fields, the 
system is started at saturation (all spins aligned) and the 
field is then oscillated back and forth while being slowly 
reduced in magnitude. For an exact demagnetization, 
the field reversal points are sufficiently reduced so that 
the last avalanche of spins does not align with the field. 
One may also perform a non-exact demagnetization by 
reducing the field by a fixed amount AH each oscillation. 

Such a field history produces an infinite set of nested 
subloops. The demagnetization curve itself is then com- 
posed of the last avalanche of each subloop, which brings 
the field back to its return point. By examining the scal- 
ing behavior of these subloops one should then be able 
to determine the corresponding behavior for the demag- 
netization curve. 

The scaling behavior of subloops has been shown to be 
dependent on that of the saturation loop plus one new 
exponent [||, ||] . The scaling forms for the saturation loop 
and subloops may be combined to give a general form de- 
pendent on the field H, disorder R, subloop return point 
magnetization M max , and the parameter of interest. In 
particular, near the critical point the avalanche size dis- 
tribution scales as 

D(S, H, R, M max ) ~ S- T V{S af5S ti, S a r, S°'e), (2) 

the correlation function as 

G(x, H, R, M max ) ~ x-W-^Gix/tir, h', e)) (3) 

where £(r,h',e) ~ r~ v y(h! /r@ s ,r~ v / v e), and the num- 
ber of spanning avalanches scales as 

N(L, H, R, M max ) ~ L°+W v M{L x l v r, L^'e, h'/r" 5 ). 

(4) 

Here h! = h + Br + Bhe is the tilted scaling axis, where 
B and B^ are non-universal constants, with h = H — 
H*, r = (R- R c )/R, and e = {M maXa - M max )/M max 
H |^, [l6| [TtJ . Note that as the last avalanche in a subloop 



occurs at M max , M max also lies on the demagnetization 
curve so that M maXc = Mf}. The critical point of the 
demagnetization curve is at (H^,M d ). Corrections to r 
and e are not significant for the following calculations. 

For e > subloops are steepest at their endpoints, 
where h! = 0. From this one may construct the 
scaling forms for the demagnetization curve by real- 
izing that at b! = the forms for subloops and 
the demagnetization curve must be identical. As an 
example, the avalanche size distribution would give 
D d (S,H,R) = D(S,H,R,M max ) at h' = or equiv- 
alently S- T *V d {S< r *P* s *h' d ,S< T *r) = S- T V(Q, S a r, S a 'e) 
where subscript d denotes quantities associated with the 
demagnetization curve and h' d — h + B d r. Notice that 
h' = implies e = —(h + Br)/Bh so that both sides have 
the same functional dependence. From this we may read 
off that one expects r d = r, B d = B, a d /3 d S d = a' , and 
a d = a. 

The integrated avalanche size distribution for the de- 
magnetization curve is obtained by integrating over the 
field h' d which is equivalent to integrating over e at h! = 0. 
Performing this integration yields 

D d nt {S,R) - S~ (Td+adt)d5d) V d nt {S' Jd r) - S- {TW) V lnt {S a r). 

(5) 

From this one obtains another scaling relation, r d + 
cr d (i d 8 d = r + a'. 

A similar analysis of the correlation function yields the 
relations v d = v, pd/v<i = {(3-/35 + a 1 / a) / v, and r] d = r/. 
Also, integrating the number of spanning avalanches over 
e yields the additional relation 6 d = 9 + ([35 — a' /a) jv. 
The associated scaling forms are 

G d nt (x,R) - x- {d+Pd '^g d nt {xr Vd ) (6) 

and 

NjP*(L, R) - L 6d N l d nt {L l/vd r) (7) 

respectively. It is possible to determine other relations, 
but from these we have already determined all the de- 
magnetization exponents from those associated with the 
saturation loop and subloops. A summary of these rela- 
tions is given in Table |. Note that all demagnetization 
quantities may be expressed as the corresponding satu- 
ration loop value plus a function of the saturation loop 
exponents times a correction factor X' = ((35— a' / cr)/v — 
(35/v- (3 d 6 d /v d (see Table |). 

The integrated avalanche size distribution was mea- 
sured on the demagnetization curve in systems with 100 3 
spins for several disorders above the critical disorder. 
The resulting distributions are shown in Figure |2|. As 
with the saturation loop, the cutoff size is found to in- 
crease as the disorder approaches the critical disorder. 
The avalanche distributions were collapsed using the scal- 
ing ansatz given by Eq. and are shown in the inset of 
Fig. S. Within errorbars, the resulting scaling function 
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TABLE I: Relations between exponents describing the sat- 
uration loop (unprimed), subloops (primed), and demagneti- 
zation curve (subscript d) scaling behaviors. 



TABLE II: Universal critical exponents found from scaling 
collapses for both the demagnetization curve and saturation 
hysteresis loop. 



Relation 


Source 


Exponent 


Demagnetization 


Saturation" 


Td ± o~dPd5d — t + a(55 — avX 1 


D in t (S, R) 


t ± apS 


2.10 ±0.05 


2.03 ± 0.03 


Pd/vd = P/v - X' 


G int (x,R) 


l/o 


3.9 ±0.4 


4.2 ±0.3 


= a 


D lnt (S, R) 


d + P/v 


3.1 ±0.2 


3.07 ±0.30 


6d = + X' 


N(L,R) 


1/v 


0.71 ±0.10 


0.71 ±0.09 


T d = T 


D(S,R,H) 





0.01 ±0.01 


0.015 ±0.015 


Vd = V 


Gint(x,R) 


pjv 


0.03 ± 0.02 


0.025 ±0.020 



X 1 = (J3S - o'/o)/u = a (pS - p d 6 d )/v 



"From relations 1, 3, and 5. 



"Reference 
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FIG. 2: Integrated avalanche size distribution curves for 
100 3 systems at disorders R = 2.5 J, 2.7J, 2.9J, and 3.3J and 
averaged over 10, 5, 4, and 1 random seeds respectively. The 
scaling collapse is shown in the inset where r — (R — R c )/R 
with R c — 2.12J. The critical exponents are (t + aP5)d = 2.1 
and l/<Jd = 4.1. 
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FIG. 3: Size of the magnetic discontinuity AM™ 1 due to sys- 
tem spanning avalanches for system sizes of 20 3 , 40 3 , 70 3 , and 
100 3 averaged over 400, 60, 80, and 10 random configurations 
respectively. Curves are collapsed using r = (R — R c ) / R with 
R c — 2.12 J. The inset shows the collapse using the critical 
exponents Pjv = 0.03 and 1/v = 0.71. 



is identical in form to that found for the saturation loop 
multiplied by a constant. The resulting exponents are 
listed in Table |j| and are consistent with the relations 
(Table |l|) predicted from the previous analysis. Spin-flip 
correlation functions also collapse according to Eq. || Q . 
The measured exponents are listed in Table Q. 

Although fid/vd was measured in the correlation func- 
tion collapse, due to its small magnitude the value cannot 
be reliably extracted from the dimension of the system 
in the equation d + Pd/^d- However, one may directly 
measure Pd/^d from AM™', the change in magnetiza- 
tion due to all system spanning avalanches p6| . For the 
demagnetization curve, the scaling form is 

AM^ nt - L-W^AMffl'^r). (8) 

The jump in magnetization was measured for system sizes 
of 20 3 , 40 3 , 70 3 , and 100 3 for a range of disorders. The re- 
sulting curves are shown in Fig. □ with a scaling collapse 



using Eq. |S| in the corresponding inset. The exponents 
are listed in Table || 

The number of system spanning avalanches in the de- 
magnetization curve was also measured for the same sys- 
tems as used in the size of spanning avalanches collapse. 
The resulting exponents are listed in Table |l[ Due to 
fluctuations the curves do not collapse well near their 
peak values. The exponent 9d, which scales the peak 
value, has thus been left with a large uncertainty and 
could possibly vanish III . 

Comparing the exponents from the two finite size scal- 
ing collapses shows that the values of 1/v are in excel- 
lent agreement with the prediction that they should be 
equal in the saturation loop and demagnetization curve. 
Similarly, the two exponents governing the power-law be- 
havior at the critical point are found to differ from the 
saturation loop values, as predicted. However the error 
bars have significant overlap. A simpler relation may be 
obtained by adding the equations for 9d and Pd/vd giving 
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FIG. 4: The demagnetization curve (solid line) and several 
subloops (dashed lines) at a disorder of R = 1.8 J and Ji n f = 
0.25J in a 100 3 system. 



from those of the saturation loop and subloops. Ex- 
periments to confirm such behavior have not yet been 
performed, but hard magnets or systems where disor- 
der dominates are expected to exhibit such behavior [ fl9| . 
Soft magnets are not expected to exhibit such behavior 
as realized when long range demagnetizing fields are in- 
cluded in the RFIM. The system is then self organized 
critical, with power law distributions of avalanches for 
both the saturation loop and subloops. 
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man, for very useful discussions. This work was sup- 
ported by NSF grants DMR 99-76550 and DMR 00- 
72783. K. D. also gratefully acknowledges support from 
an A. P. Sloan fellowship. We also thank IBM for a gen- 
erous equipment award and the Materials Computation 
Center for the use of their IBM Workstation Cluster. 



0d + fid/vd = + (3/v which is satisfied by the measured 
exponents. The value of the correction factor may be 
estimated from the finite size scaling exponents giving 
X 1 = —0.005 ± 0.005. The power law exponent from the 
avalanche size distribution also contains the correction 
factor, but due to relatively small system sizes, finite size 
effects do not allow one to reliably extract X 1 . While X 1 
was found to be slightly negative, because of the large 
uncertainty in the critical exponents for both the satura- 
tion loop and demagnetization curve, one cannot rule out 
the possibility that X 1 vanishes (as occurs on the Bethe 
lattice [^0) ) . If such were the case then all the critical ex- 
ponents would be identical for both the saturation loop 
and demagnetization curve. 

In the above analysis we have ignored long-range de- 
magnetizing fields. To examine the effect of such fields 
one now sets Ji n f > in the Hamiltonian, Eq. |l|. A de- 
magnetization process may again be applied to the sys- 
tem by applying an oscillating field and decreasing it by 
a given AH each oscillation. The resulting curves for a 
100 3 system at a disorder of R — 1.8 J with J in f — 0.25 J 
are shown in Fig. For clarity only a few subloops 
are shown. All the curves, including the demagnetiza- 
tion curve, show a linear M{H) behavior as the field is 
increased, as one would expect for a system undergo- 
ing domain wall propagation. Additionally, for the same 
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law distribution of avalanche 



cur\es one nnas a power law custriDution 
sizes, as expected for a self organized critical system 
H |(|. Preliminary experimental results by A. Mills and 
M. B. Weissman using a soft magnetic material suggest 
the results of Fig. || are correct |l8) . 

In the absence of long range demagnetizing fields, the 
demagnetization curve of the RFIM exhibits disorder in- 
duced critical scaling in quantities associated with its 
Barkhauscn noise. Additionally there are no new in- 
dependent critical exponents, rather all may be derived 
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